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ABSTRACT 

Forthcoming astrometric missions such as the Gaia satelHte will bring to the fore the problem 
of estimating the enclosed mass from a set of positions, radial velocities and proper motions of 
tracer stars. Here, we show how to construct the tracer mass estimator when the proper motion 
data are available in addition to the usual line-of-sight velocity data. Notably, the mass esti- 
mators do not require any assumption on the anisotropy, as it is naturally incorporated through 
the different components of proper motions. In addition, the separate treatment of the proper 
motions and the line-of-sight velocities is desirable because they are observationally indepen- 
dent and thus the propagation of the combined uncertainties is rather straightforward. The 
extension to projected data is also sketched, together with a possible application of measuring 
the masses of Galactic globular clusters. 

Key words: galaxies: fundamental parameters - galaxies: general - galaxies: haloes - galax- 
ies: kinematics and dynamics - dark matter 



1 INTRODUCTION 

The great era of Galactic Cartography is upon us. The last twenty 
years have seen ambitious photometric and spectroscopic surveys 
of the Galaxy provide magnitudes and radial velocities for many 
millions of stars. Astrometric surveys of the Galaxy began with 
the n ow comparatively modest Hipparcos mission (see IPerrvman I 
I2OO9I for a summary of its achievements) and will culminate with 
the launch of the Gaia satellite in 2013. This will provide micro- 
arcsecond astrometry for ~ 10^ stars brighter than G ~ 20 
JPerrvman et alj|200ll : iTuron et"al] |2005'), as well as radial veloc- 
ities for stars brighter than G ~ 17. Positional and kinematic data 
on bright stars in the Galaxy and its neighbours will be routinely 
available by the end of its mission. 

This wealth of new data motivates the study of new theoretical 
questions. A classical problem in dynamical astronomy is to esti- 
mate the enclosed mass from the relative positions and velocities 
of a collection of bound test particles or tracers. Hitherto, work 
has primarily concentrated on the case in which the velocities are 
line-of-sight or projected. Here there has been an extensive body of 
work starting with methods essentially e xploiting the virial theorem 
jPagel I1952I : iLimber & Mathewd 119601). The problems with suc h 
virial estimators were enunciated by Bahcall & Tremain3 ( Il98lh . 
who introduced the notion of projected mass estimators based on 
sums of the projected mass contributions, namely [(projected dis- 
tance) X (radial velocity)' / G]. In this form, the estimator is natu- 
rally suited to the case of a point mass (the Kepler potential). The 
extension to scale-free potentials with ijj oc r^" was developed by 



I White! ( Il98lh . lEvans et al.l j2003l) and IWatkins et all ( I2OIOI here- 
after WEA). The estimator then takes a different form, namely sums 
of cont ributions of the form [(distance)" x (radial velocity)^ / G]. 
Finally, lAn & EvansI j201 iL hereafter Paper I) have shown how to 
adapt these ideas for an arbitrary spherical potential by choosing a 
weight function. 

The case of proper motions has so far received very little atten- 
tion. The coming of age of Gaia will make the availability of proper 
motions and radial velocities standard for all objects in the nearby 
Universe brighter than G « 20. Hence, this paper is devoted to the 
problem of estimating the enclosed mass from a tracer population 
when their space velocities are known. 



2 TRACER MASS ESTIMATORS 

Throughout this paper, we operate under the simplest assumption 
that the gravitational potential is spherical. A tracer population of 
stars with a number density v resides in the potential. This is sam- 
pled to provide A'to, objects with known positions, radial velocities 
and proper motions. We shall usually consider the application to 
the problem of estimating the dark halo mass within the radius of 
the outermost datapoint TaM- 
We start with the spherical leans equation, which reads 



d(vo-^) 
dr 



where 
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IS the lBinnevldigS^ anis otropy parameter. This can be integrated 
into the form (see Paper j eq. 11) 



GM„„, = {[3 + a(r) - 2f}(r)] -^) - 3r„„r, 



(3) 



where the angled bracket represents the average over the whole 
tracer population. Here Mout is the total halo mass within radius 
J-out, and 



M(r) 



(0 < r < r„„,) 



(4) 



is the halo mass profile M(r) normalized by the total mass such that 
/i(''out) = 1> whereas 



dlog[r 'M(r)] dlog/i(r) 
a{r) = 7-: = 1 - 



dlog r 



dlog r 



(5) 



is the effective index for the gravitational potential. Finally, the sur- 
face term is given by 



v(rou,)cr^(rou,) 



3M„ 



> 'out — . T • 



(6) 



Note that Vou, is the mean number density of the tracers in the sphere 
of the radius of ra^t- 
Next, we observe from equation ^ that 



\^) = ir ^dr= 2(1-/?)^ 
\ I NtotJo 1^ \ ti 



(7) 



where = (Ug + o^)? is the tangential velocity and tr, = (a^ + 

cr^)^ is the corresponding velocity dispersion. Hence, equation l[3j 
is equivalent to 



lrv:\ , 
GM„„, = / + ( y ) - 



(8) 



which is our desired result expressing the mass Mom in terms of 
averages of observables. 

The same can also be derived by directly integrating equation (T} 
whose left side reduces to 



d(vo-J) 



1 d(i^va-f) va\ 



(9a) 



dr r or r 

With an arbitrary weight h = h{r), this then results in 

f™' ,d(rVcrJ) r»' , , f""' , 

Arhr— drhvo-^r^ = - drhvGMr^. (9b) 

Jo dr Jo Jo 

The first integral here reduces to 



Jo 



dr hr 



2d(rWJ) 
dr 



' 2 ^ dihr') 

drr vaz 

' dr 



Jr-rc 


hrva--]^ - dr hvcrtr (9c) 

'0 Jo dlogr 



through integration by parts. Given that the inner boundary term 
hr'*vcr^\Q = vanishes, replacing the first integral in equation i9b\ 
by equation l [9cl ) and multiplying the result by 4nNfJ^ , we find that 



(hGM) 



We recover equation ([8]l by choosing h = ft ' since d log h/d log r = 
-dlog///dlogr = a - 1. 



2.1 Power-Law Weight 

2.1.1 Scale-free potentials 

lEvans et"ail ( l201ll . hereafter BAD) have demonstrated that for most 
practical purposes adopting the scale-free assumption is a good ap- 
proximation for choosing the proper mass estimator. This amounts 
to assuming that the mass model is given by 



ii{r) = (— )' " (0 

\ ''out / 



< f < '■out), 



(11) 



and so a = a is constant. The exact models of the Kepler poten- 
tial (a central point mass), the logarithmic potential (a truncated 
singular isothermal sphere), and the harmonic potential (a uniform 
homogeneous sphere) correspond to the cases that a = I, a = 0, 
and a = -2, respectively. With the mass model of equation illl . 
equation l[8j reduces to 



= (l+a)< 



(12) 



Given that the tracer density profile essentially vanishes at finite 
Toui, the surface term should drop out (unless they ar e press ure- 
confined with a hard cut-off). Alternatively, following IwEAI the 
surface term may be explicitly solved under the assumption of ex- 
act power-law behaviour for the tracer density (v oc r"'') and the ve- 
locity dispersions, although strictly speaking the resulting solutions 
cannot be physical everywhere. It follows that yS is constant since 
the Jeans equation implies that crj, c^^, and GMr^^ all follow the 
same power-law. Therefore we can solve for v(rou,)/Vou, = |(3 - y) 
and crj = (y-2fH+ay'GMr^\ If we substitute these for the surface 
term, then we find that: 



+ 3f^ 



3-2/3 + 0- GMo, 



'■out y-2p + a r„ut 

Together with (vfr") = 2(1 - /3){v^.r"), equation | |12| | results in 
GM, 



(13) 



(a + 7 - 2) (vy) + Ur") 



(14) 



This is equivalent to equations (16) and (24) of IwEAI if the 
anisotropy parameter is given by /3 = 1 - {vy)/(2{vjr°^)). 



2.1.2 General cases 

In reality, the scale-free mass model is a rather artificial abstraction. 
Never the less, equation JlOt is actually valid for any mass profile 



and arbitrary weights. Hence with h 



and further assuming 



that the surface term vanishes as rout ^ oo, we have that 



G{r"-'M) = (1 + a)U.r") + Ur" 



(15) 



We have used the symbol a to emphasize the fact that this needs 
not be the same as the power-law index of the potential. That is to 
say, the scale-free estimator on the right-hand side actually returns 
the radial moments of any given (not necessarily scale-free) mass 
profile. Although it is clear that the left-hand side is directly related 
to the total mass, the specific constant factor associated with it in 
general depends on both the mass profile and the tracer density 
profile. There are some exceptions that are worthy of note here. The 
first is the case of the power-law mass profile with a = a, for which 
(r""'M(r)) = <r|j'7,'Mout) = Cu'^out is constant regardless of the 
tracer density profile. The second is the self-consistent gravitating 
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Table 1. Mean values {vlr") for satellites in Table 1 of lWEAl 



''out 


300 kpc 


200 kpc 


100 kpc 






(xlO'i M©) 






2.86 ±0.87 


1.75 ±0.40 


1.20 ±0.22 


(w/o Leo I) 


2.09 ± 0.42 






(w/o Leo 1, Her) 


1.76 ±0.27 


1.45 ±0.28 





Here a = 0.55 following I WEAI Only satellites within r < Tom are counted. 
These consist of 14 satellites with r^^t = 100 kpc, 21 (including Her) with 
rout = 200 kpc, and 24 (including Leo 1 and Her) with rout = 300 kpc. The 
quoted uncertainties, which range from ~ 15 to ~ 30%, are purely due to 
the scatter, i.e., cr/W''^ = sj yjN — 1 where cr^ and s- are the unbiased and 
biased sample variance estimate and N is the number of data points. 



system for which the a = \ case results in 

{M)= drv?M= dr 

Mot Jo A^toi Jo d'- 



M 



2Mto,Jo 



dM^ 
dr 



Moi 
2 



(16) 



because dM/dr = 47tr-p where p is the mass density and the ratio 
p/v = Mtoi/A'to( is constant. (By contrast, for a central M. point- 
mass, we would have (M) = (M.) = M,.) An additional general 
result deduced from equation l ll5t is for 5 = 0, that is, 

{^)-{^).{^) = i.\ (17) 
which is simply the scalar virial theorem for a spherical system. 



3 THE MASS OF THE MILKY WAY 

3.1 Present Proper Motions of the Satellites 

IWEAl have used combinations of the averages in the form of {v^^r") 
and (u^r") = (3 - 2[}){v^.r") to estimate the mass of the Galaxy 
from the kinematics of its satellites including their proper motions. 
Whilst their approach of incorporating proper motions has a valid 
justification in that the anisotropy value derived by combining the 
line-of-sight (los) and proper motions appears to be severely biased 
towards tangential motions, it is still instructive to use the consis- 
tent mass estimator without a priori assumptions on the anisotropy. 

Here we use equa tion (I14t to combine the los and proper mo- 
tion data gathered in lWEAl and estimate the Galactic mass. We first 
note the scaling relation for observable combinations of kinematic 
properties ((las a"' = micro-arcseconds per year): 



Gr\ 



out 

a 



2.325 • 10 



5.225 • 10' 



'■out y~'Y r \"l ve \^ 
100 kpc / 1 100 kpc M 10 km S-' / ' 



n , / ''out 

^^Took^ 



I 1 100 kpc/ \10^asa-'/ 
(18) 



Here, vc and ii are the los velocity and the proper motion of the 
satellite, respectively. For the satellites considered here, the Galac- 
tocentic distances are large enough to ignore the sun being off- 
centre in the Galaxy. We shall therefore assume that Vy = Vf and 
D, = r/x, and the use of equation J14t is appropriate. 

Table[T]lists the contribution to equation il4\ of the radial motions 
of satellites assembled in Table 1 of WEA with a = 0.55, which 
is the 'best-fitting' value of the power-law index to the potential 



Table 2. Numerical contribution to the pro per motion mass estimator 
G^'r^-f^u^r" for satellites in Table 4 of lWEAi 



''out 


200 kpc 


100 kpc 




(xlO" 


Mo) 


Car 


5.3 ± 3.6 




Dra 


90.6 ±47.1 


66.3 ± 34.5 


For 


60.6 ± 9.4 




MCs 


48.5 ± 2.4 


35.5 ± 1.8 


Scl 


0.43 ± 1.20 


0.31 ±0.88 


Sex 


4.1 ± 12.2 


3.0 ± 9.0 


UMi 


10.9 ± 6.7 


8.0 ± 4.9 



Here again a = 0.55. All listed satellites have r < 100 kpc except for Ca- 
rina (r = 102 kpc) and Fornax (r = 140 kpc). For the Magellanic Clouds, 
we choose r = 50 kpc for the barycentre (assuming the mass ratio of 
LMC/SMC ~ 11). The quoted errors are the directly propagated uncertain- 
ties of the proper motion measurements assuming no correlation between 
the two components of the proper motions. 



of a typical NEW halo as determined bv IWEaI The uncertainties 
quoted are purely statistical and typically of about ~ 20 per cent, 
comparable to ~ 25 per cent derived via Monte Carlo simulations in 
|WEA. Assuming isotropy and adopting the tracer power-law slope 
7 = 2.6, the Galactic mass based on these data can be estimated by 
multiplying these by a factor of 3. 15, t he resu lts of which basically 
reproduce those reported in Table 3 o f lWEAi 

Table |2] on the other hand provides us with the individual con- 
tributions to the mass estimator of the proper mo tions m easured 
for seven Galactic satellites reported in Table 4 of lWEAi In con- 
trast to the los velocities, the uncertainties associated with each 
measurement are rather large. This not only contributes additional 
uncertainties in the final mass estimate, but also complicates the 
proper weighting scheme appropriate for the average related to the 
mass estimator. In Table [3] we present the result calculated with 
equation jl4t . Eor the average of the proper motion contributions, 
the results based on both the unweighted arithmetic mean and the 
mean weighted by the reciprocal of the observational variance (i.e., 
cr"' where cr is the associated uncertainty reported in Table |2j are 
given. Compared to the results based on the los motions alone, the 
final statistical uncertainties are quite larger ranging from ~ 50 per 
cent to a factor of two. More seriously, the results vary in signif- 
icant amounts depending on the weighting scheme and the choice 
of outliers (which may be due to observational errors and/or true 
unbound objects), although they are formally consistent with one 
another given large error bars. It appears that the safest conclusion 
that can be drawn at this point is that the existing proper motion 
data do not provide any meaningful discriminating information on 
the mass of the Milky Way without any a priori assumption on the 
motions of the Galactic satellites as a whole or the particular target. 



3.2 Tlie Age of Gaia 

Let us now estimate the benefits that Gaia will bring. We do this by 
constructing mock datasets of los velocities and proper motions for 
sets of dwarf galaxy satellites extracted from a suite of simulations, 
the Galaxies-Intergalacti c Medium Interact ion Calculation, which 
are described in detail in ICrain et alj ( |2009|) . They consist of a set 
of hydrodynamical resimulations of five nearly spherical regions 
(~ 2 0/i~' Mpc in radius) extracted from the M illennium Simula- 
tion jSpringel et alj|2005h . loeason et al.l ilOl ll) extracted a set of 
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Table 3. Mass of the Milky Way 



'"out 






200 kpc 

(w/o Her) 


100 kpc 


unweighted 


33.5 


± 


15.4 


33.2 ± 


15.4 


24.0 


± 


14.9 


(w/o Dra) 


23.7 


± 


11.1 


23.3 ± 


11.1 


13.1 


± 


8.6 


(w/o MCs) 


30.7 


± 


17.9 


30.4 ± 


17.9 


20.8 


± 


18.8 


(w/o Dra,MCs) 


18.3 


± 


11.9 


18.0 ± 


11.9 


5.2 


± 


4.8 


o""- -weighted 


12.1 


± 


10.6 


11.8± 


10.6 


8.6 


± 


10.4 


(w/o Dra) 


12.1 


± 


11.4 


11.8± 


11.4 


8.6 


± 


11.6 


(w/o MCs) 


4.1 


± 


6.5 


3.8 ± 


6.5 


2.0 


± 


4.9 


(w/o Dra,MCs) 


4.1 


± 


6.8 


3.7 ± 


6.8 


2.0 


± 


4.2 



The estimator uses equation )14t with a = 0.55 and y = 2.6. The mass is in 
the units of lO" Mq. The uncertainties are purely statistical incorporating 
both scatters and the measurement uncertainties found in Table|2] 



galaxies that resemble the Milky Way. The catalogue consists of 
431 parent haloes and 4864 associated satellite galaxies. 

The mass estimators provide the total mass within the radius of 
the farthest tracer (« r^M)- We compute the 'true' mass within the 
virial radius of each halo and compare to the masses found via our 
estimator. We use all satellites, but check that our results are not 
significantly affected when only luminous satellites are included. 
Figure [T] refers to data sets of radial velocities and true distances 
(dashed histogram), and radial velocities, proper motions and pro- 
jected distances (solid histogram). An uncertainty of 100 H-as a"' for 
each star is assumed; for comparison, Gaia's tar geted accuracies 
are 10 jxas a^'at G ~ 15 and 150 (xas a^'at G ~ 20 jLindegren et al.l 
l2008h . For each tangential velocity measurement, a Gaussian dis- 
tribution centred on the 'real measurement' with error is derived 
from the proper motion (and so depends on distance). For each 
satellite, it is assumed that Gaia will detect 100 stars, which is 
reasonable enough for the case of the satellites of the Milky Way 
jWilkinson & Evansll 19991) . The left panel of Figure[T]uses the es- 
timator given by equation l ll4t . the right panel uses equation l ll2t 
neglecting the surface term. 

We quantify the improvement by m eans of two statistical mea- 
sures as originally introduced bv lEADt First, we define the Fraction 
of Reasonable Estimates (FRE) as the fractio n of estimates withi n 
the factor of two of the true mass (see also Deason et al. I l201lh . 
We also give the Inter Quartile Range (IQR) of the mass estimates, 
which gives a good indication of the spread. We see that at Gaias 
accuracy, the proper motions can contribute significantly to an im- 
proved estimate of the mass of the Milky Way. Even with the small 
number of satellites, the inclusion of proper motion data via the 
estimator of equation J14l > means that 93 per cent of estimates are 
within a factor of two, as opposed to 82 per cent with radial veloci- 
ties alone. The solid histogram is also more narrow than the dashed, 
as indicated by a smaller IQR of -0.05 < logiM^^t/Mt^u^) < 0.09. 
Notice from the right panel that the estimator of equation il2\ also 
performs very well with regard to the FRE, but the neglect of the 
boundary term c auses an offset so the estimator is biased. 

As is noted by'EAD', there are a number of ways in which mock 
catalogues from simulations differ from the assumptions used to 
derive the estimators. For example, dark haloes are not generally 
spherical, infall continues to the present day, and the observed satel- 
lites are not necessarily virialized and well-described by an equilib- 
rium distribution. Therefore the performance of the mass estimators 
is very encouraging. 



4 PROJECTED DATA WITH PROPER MOTIONS 

In most practical cases, the kinematic data available are limited by 
projection. Hence, the naturally observable velocity dispersions are 
those along the line of sight and of proper motions. In addition, for 
most dynamical systems of interest, their spatial extents are typi- 
cally much smaller than their distances from us as a whole. On the 
positive side, this indicates that, from the observer's perspective, 
the observed los velocities V[ basically correspond to the coordi- 
nate velocity component v- if we have chosen the cylindrical polar 
coordinate centred at the centre of the tracer population and the pos- 
itive z-axis running from the centre to the observer. In addition, the 
remaining two orthogonal components of the velocity Dp = (vr, v^) 
are linearly proportional to the two components of the proper mo- 
tion with v^, = sfi where s is the distance to the system as a whole 
and fi is the proper motion vector. 

On the negative side however, the radial distance of any individ- 
ual tracer to the centre of the system is very difficult to determine 
precisely. In particular, since the absolute uncertainty (Sr) of the 
radial distance (r) measurement is in the same order as the uncer- 
tainty (Ss) of the individyal distance (s) measurement, i.e., 6r ~ 6s 
(NB., the eiTor in the distance measurement typically dominates 
that in the angular separation whilst the precise determination of 
the distance to the centre only improves the precision by ~ V2), the 
fractional uncertainty of the radial distance is roughly estimated to 
be (6r/r) ~ (6s/s)(s/r) ~ (Ss/s)/ip where (p ~ r/s is of order 
the angular extent of the system. In other words, ~ 10 per cent 
estimates of the radial distances for tracers in a degree-sized sys- 
tem require the individual distance measurement in a better than 
~ 0.2 per cent precision. Consequently, the distance handle associ- 
ated with the tracer is likely to be limited to the projected (linear 
or angular, depending on whether the distance to the system as a 
whole is known or not) separation. In the follo wing, we develop the 
projec ted rt iass estimators simi lar to those of ISahcall & Tremaind 
l ll98lh and lHeisleret"ai] ( Il985l) using V( and R as observables but 
further augmented by including the transverse velocity »j, derived 
from the pro per motion. We note that this problem has been partly 
addressed bv lLeonard & Merritt H98a) . 



4.1 The Projected Mass Estimator with Transverse Velocities 

Let us first consider the coordinate system with its origin located 
at the centre of the system and the positive z-axis running from the 
origin to us. Then, the tracer velocity components in the spheri- 
cal polar coordinate (u, , vg, v^) transform to those in the cylindrical 
coordinate (%, v^, v-) according to 

Vf: = VrSind + Vf)Cos9; v- = VrCosd - vgsinO. (19) 

If the tracers are described by a generic two-integral distribution 
that is non-rotating and spherically symmetric, then the velocity 
dispersion tensor is diagonalized in the spherical polar coordinate 
with CTgg = a^^. These transform to the velocity dispersion tensor 
in the cylindrical polar coordinate via 

crl^ = a-j sin^e + crj cos^e = (1 -/3cos-e) a] 

al, = al = {l-p)a] ^^^^ 

al, = a], cos-e + cr-g sin^O = (1 - /? sin^O) cr\ 

a\, = (crj - o"^) sin 6 cos 6* = fS sin 9 cos 6 cr'^, 

where = crj^ and crj = a^g, and yS is given by equation l|2j. 
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Figure 1. Left: Histograms showing the ratio of the estimated mass to the true mass for the two different mass estimators. The dashed histogram uses sets of 
radial velocities and distances with the tracer mass estimator of lWEAl The solid histogram uses sets full space velocities and distances, with mass estimator of 
eq. U4) . The FRE and IQR are the fraction of reasonable estimates and the inter-quartile range as defined in the text. The means of the distributions are shown 
by vertical dotted lines. Right: As the left panel, but now the solid histogram uses eq. )12t neglecting the surface term. 



Using equation l |20| l, it is trivial to find for a > -2 that 
vlR") = drrv(r) dd sin cr^^r" sin"0 

Mot Jo Jo 

= B,,[(2 + aXv'X) + (v'gr''}] 
vlR") = (3 -H a)BMr") 
vlR'') = B„[{vy) + {2 + a){vlr'')\ 



(21a) 

(21b) 
(21c) 



where B^ ;, = T{a)T(b)IT{a + b) and T{x) are the beta and gamma 
function and 



V^r(i + f ) 
4r(^) 



(21d) 



The result extends to the a = -2 case (homogeneous sphere) with 



{4&(R')) = dz j ARRyal^b(R^ 



(22a) 



dz vtr' 



^^Ir=o.c 



where 8(x) is the Dirac delta. Since {R = Q,z > Q) and {R = Q,z < 
0) correspond to (r = z, f = 0) and (r = \z\,9 = n), we find from 
equation i20\ that 



2n 



Jo 



(22b) 



{.lb(R-))=^rdr. 

'Vtot Jo 

{vlbiR^))=^rdrya^^ 

'Vtot Jo 

This result is identified with the formal limit of equation J2U as 
a -» -2+ once we note that lim„^_2+ R" IT{\ -i- f ) = HR-)- 

For an observer located sufficiently far away, the axial and sec- 
tional velocities are basically the los and the transverse velocities. 



respectively. It follows that (v^R"^ = (v^R"^ and 

{vl,R") = B„ [(2 + aXv^j-") + (4 + aXv^r"}] . (23) 

If the two components of the transverse velocities are separately 
measurable, we can derive the right-hand of equation l |15t from 

(l+a){vlR") + {vlR'') = (2+a)B,[(l + «)(.>") + 2<.|r")] . (24a) 

On the other hand, the same may be found by combining the los 
and the two-dimensional transverse velocity dispersions 

{i/^,R") - {v^R") = B„ [(1 + a){vy) + 2{vlr'')\ . (24b) 

Directly substituting these in equation l |15l l then leads to 



{l+a){vlR") + {vlR") 
2 + a 



(vlR^yUR'A, (25a) 



which is again valid for a > -2. We note that equation l l25at with 
a = -1,0, 1,2 are already found in .Leonard & Merritt (1989 ). If 
a = -2, then strictly the average is not convergent, but using equa- 
tion l l22t . we find the proper limiting equation, namely. 



{G{r-'M) = {vlb{R^))-Ub{R^)) 



(25b) 



whereas (vj^b(R^)j = (v^b(R-)y Provided that the distance s to the 
system as a whole is known independently, equation J25at is the 
proper extension of the projected mass estimator to include the in- 
formation from proper motions. Note that, if the vectorial disper- 
sion of the proper motions is known, incorporating the los velocity 
dispersions in fact over-constrains the system under the assumption 
of the two-integral spherically-symmetric distribution. 

4.1.1 Mass measurements of globular clusters 

We can test these formulae in the context of applications. A pos- 
sible class of targets would be nearby dwar f sphe roidal/elliptical 
galaxies, similar to the examples found in IEADI However, the 
proper motion dispersion for an imaginary satellite of mass 10' Mq 
and radius 0.5 kpc located at 50 kpc is about ~ 40 ^.as a"'. Even 
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Figure 2. Left: Estimates of enclosed mass of a globular cluster using 10 (dashed) and 100 (solid) tracers from lO' Monte Carlo realizations. Proper motion 
errors of 100 iJ.as a" 'and a distance of 5 kpc are assumed. Right: Estimates of enclosed mass as a function of globular cluster distance, namely 5 (solid), 20 
(dashed) or 50 (dotted) kpc, assuming 100 stars can be measured per globular cluster. Model has parameters a = and y = 3. 



with Gaia, these will be challenging observations even if there ex- 
ist large enough samples of sufficiently luminous tracers, which 
is itself unclear in any case. Alternatively, we can consider mass 
estimates of Galactic globular clusters. In fact, for a few bright 
clusters, t here already exist very good inter nal proper motion data 
sets (e.g.. [Anderson & van der Marelll2010 '). Yet the existing data 
for these are sufficiently large and of high quality so that the 
complete dynamical modelling is technically feasible (see e.g., 
Ivan der Marel & Andersonl2010l) - note that there is no degeneracy 
involving anisotropy with the proper motion data incorporated. In 
principle, we may still utilize these for the projected mass estimator 
but its application is clearly didactic at best. Instead we investigate 
Gain's capability regarding possible applications of projected mass 
estimators for 'typical' globular clusters in varying distances. 

In contrast to dwarf satellites, the main challenge of globular clus- 
ter observations is expected to be crowding, and so the choice of 
the tracer subsample is important for any practical application of 
mass estimators. Blue horizontal branch (BHB) stars have a G-band 
magnitude of about ~ 19 at 50 kpc. So, Gaia will be able to detect 
a BHB star out to 50 kpc and measure their proper motions with 
accuracies of better than 100 (xas a"' . By contrast, radial velocities 
will usually be missing, as Gaia's radial velocity spectrograph has 
a limiting magnitude of G « 17. Most g lobular clusters ha rbour 
significant populations of BHB stars (e.g., lBrown et al.ll2005l) . and 
so equation l l25at can be exploited to measure the mass of globular 
clusters in the Milky Way from proper motions provided by Gaia. 
To test the feasibility, we use the the self-consistent distribut ions of 
power-law densities wit h p oc r"' in logarithmic potentials jWhitd 
1 19811 : 1 Evans et al. l fT997l) to generate mock datasets of proper mo- 
tions and projected distances, degraded with the Gaia errors. The 
estimated mass computed via equation l l25at can then be compared 
to the known true mass. Figure |2] shows the results of Monte Carlo 
simulations. In the left panel, the globular cluster is assumed to lie 
at a distance of 5 kpc and the sample size is 10 or 100 tracers with 
projected positions and proper motions. As the histograms are quite 
tight, we have replaced FRE with the Fraction of Good Estimates 
(FGE), where "Good" means that the mass estimate is within 10 per 
cent of the true mass. We see that 10 tracers are too small for re- 



liable results, but 100 tracers are ample. The right panel shows the 
estimated masses as a function of distance, with the globular cluster 
assumed to lie at 5, 20 and 50 kpc respectively. In each case, 100 
tracers are used, and we see that the results are still acceptable at 
distances of 20 kpc, with 67 per cent of the mass estimates being 
good. With increasing distance, however, the accuracy degrades so 
that by 50 kpc, only 24 per cent of the estimates are good. 

4.2 The Projected Mass Estimator with Proper Motions 

Now we consider the case that the distance s is also unknown. It 
can in principle be measured by matching the los velocity and the 
proper motion dispersions. Specifically, the last equality of equa- 
tion l l25at is rearranged to 

(2 + a) = {vlR") + (l+a) {uIr") . (26) 

With R = s(p and By = sji where if is the angular separation from 
the centre, the distance s is found from direct observables as in 



{2 + a) 




) 







(27) 



where yu^ and /i^ are the radial and azimuthal projections of the 
proper motion. Note that the a = case reduces to s-{^-) = {vf) 
irrespective of the anisotropy, which corresponds to the classical 
statistical parallax technique of measuring distance. This is because 
only in the a = case are the velocity distributions independent of 
position, which usually assumed in statistical parallax applications 
(see e.g., iTrumpler & Weaveilil953h . 
At the same time, we can also solve for the anisotropy parameter 



(3 + a)[(/.2^'') 


- 




(3 + «) 




- 





(28) 



The unique mass estimator that follows this is 



GB„{r''-'M) = 



a + a) 




+ 






+ (14- 


a) 
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with (v^fR") = s"{vj( fi''). This again is in fact equivalent to equa- 
tions (24) and (25) o f IPaper j with yS replaced by equation J28t and 
the boundary term dro pped. I f one preferred the explicitly solved 
boundary term used in IweaI one might use their equations (26) 
and (27) with p replaced by the same equation i28t . However, we 
note that estimators in equations l[27]l-l[29j are likely to be biased if 
the sample means are used to estimate for the averages in the right- 
hand sides as these formulae involve a division by a mean whose 
corresponding dispersion is finite non-zero. 



5 CONCLUSIONS 

Hitherto, proper motion data have not normally been used in mea- 
surements of the mass of nearby galaxies like the Milky Way and 
M3 1 . Proper motions are only available for a handful of satellites 
- namely, Carina, Draco, Fornax, Sculptor, Sextans, Ursa Minor 
and the Large and Small Magellanic Clouds for the Milky Way, 
together with ICIO and M33 for M31. As these are difficult mea- 
surements, the errors are often quite large. Additionally, agreement 
between diff erent investigators regarding the proper motions for 
Urs a Minor jScholz & Irwinlll994l: jPiatek et alj|2005l) and Sculp- 
tor jPiatek et al.ll2006l ; IWalker et al.l2008l) is poor. Therefore, there 
has been little attention given to their exploitation. 

This situation is about to change. The Gaia satellite, to be 
launched in 2013, will measure the positions and proper motions 
of ~ 10' stars in our Galaxy and other members of the Local 
Group. The proper motions have a targeted accuracy of between 
10 iJ.as a"'(at G ~ 15) and 150 jxas a"'(at G ~ 20). This therefore 
motivates detailed study of the use of proper motions in mass esti- 
mators. Particularly in the case of the Milky Way, the proper mo- 
tions are crucial as the line-of-sight velocity closely coincides with 
the radial velocity as measured from the Galactic Centre. Hence, 
proper motions provide new information on the velocity anisotropy 
that is not available in the line-of-sight velocities alone. 

Here we have shown how to solve the classical problem of esti- 
mating the enclosed mass from the positions and space velocities 
of a set of tracers. The main result is equation ([8}, which shows 
how to compute the enclosed mass in terms of averages over the 
observable velocity components for a given halo mass model. If, as 
is reasonable, we use a scale-free assumption for the dark halo, then 
the mass estimator reduces to equation l ll4t . Note values of o- ~ 0.5 
for the power-law index for the potential are good approximations 
to the NFW halo jPeason et alj201 ll : lEvans et al . 201 1). 

We have used these new formulae to demonstrate that the existing 
proper motion data for the Milky Way satellites provide no mean- 
ingful discriminating information on the mass of the Milky Way. 
This though will change in the era of Gaia. Using mock catalogues 
of satellites extracted from the simulations, we have shown that our 
new mass estimator, when equipped with proper motions with the 
targeted Gaia accuracies, significantly outperforms the estimator 
used with line-of-sight velocities only. We have also given the ex- 
tension of the formulae to the instances where projected positions 
and proper motions are available. These formulae provide a feasi- 
ble method for the Gaia satellite to measure the internal masses of 
Milky Way globular clusters lying within 20 kpc. 

As is usual for a spherical system, we have assumed no rotation 
in the tracer population for our derivation of the mass estimators. 
(Clearly, any rotating population cannot be spherically symmetric 
although it may be spherical in density.) However, the tracer veloci- 
ties in our formulae appear as their second moments rather than the 
dispersions/variances (in the rest frame of the potential). That is, 



once the bulk linear motion of the tracers is subtracted, the estima- 
tors actually do not distinguish between the random and coherent 
motions of the tracers. Hence, we argue that the formulae as they 
are are still valid even if the population is in fact supported through 
the mixture of the pressure and the rotation. 
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